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ABSTRACT 


A key parameter in the design trade-offs made during AlGaN/GaN HEMTs 
development is the channel temperature. An accurate determination can generally only 
be found using detailed software; However, a quick estimate is always helpful, as it 
speeds up the design cycle. This paper gives a simple technique to estimate the channel 
temperature of a generic AlGaN/GaN HEMT on SiC or Sapphire, while incorporating the 
temperature dependence of the thermal conductivity. The procedure is validated by 
comparing its predictions with the experimentally measured temperatures in devices 
presented in three recently published articles. 
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I. Introduction 


At the present time many research groups are studying AlGaN/GaN HEMTs for 
high frequency (up to Ka-band), high power applications [1], [2]. The large power 
density generated in power amplifiers causes considerable self-heating, and a reasonably 
accurate estimate of the channel temperature is often desired. Knowledge of the 
temperature is essential as both the carrier mobility and reliability suffer. The mobility 
reduces with increasing temperature as (1/T)**2.3, with resulting decrease in DC and rf 
performance [3]. Both short and long term reliability are key parameters, which also 
reflect the need to rapidly calculate the temperature profile in a particular design. 

A self-consistent solution for the thermal field in a particular device can be 
obtained using either 2-D or 3-D simulators available in the SILVACO code [4]. Some 
recent analytical articles [5], [6], [7] give techniques to quickly ascertain junction 
temperatures under certain conditions. This article demonstrates a procedure to estimate 
the channel temperature in AlGaN/GaN HEMTs on either SiC or Sapphire, using simple 
analytic expressions. The device is modeled using the thermal resistance method, while 
considering the temperature dependence of the thermal conductivity. 



II. Analysis Procedure 


A generic HEMT is shown in Figures 1 and 2. The thermal conductivities for the 
many material layers, were determined by making rough averages of the values reported 
in the literature. Figure 3 presents the simple thermal resistance stack. The heat sources 
are thin strips in the gate-drain channels. The technique consists of determining the 
thermal footprints of the heat sources at each successive lower layer, and then 
determining the thermal resistance. 

In principle, in any heat transport problem, one solves for the temperature field 
from the heat flow equation 


F'f *i(T)VT} = - p 

Where (T ) is the temperature dependent thermal conductivity in the i-th layer, T is the 

temperature, and P is the dissipated power density. The equation is altered by defining an 
equivalent artificial temperature via the Kirchhoff transformation 



where T 0 is the reference temperature, T- Q is the transformed temperature in the i-th 
layer, and T is the real temperature there. The transform linearizes eq.(l) with the result 


*i(To) = -P 
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Our first task is to determine analytical expressions for the temperature dependent 
thermal conductivities of the GaN, SiC, and Sapphire layers. We have not found an 
analytical expression for the thermal conductivity of GaN, but from [8], [9] we may 
assume it to be 



where we have assumed the behavior of GaN is similar to that of GaP. For SiC we 
choose [10] 


r SiC (T) — 3.4- 





and for Sapphire we use the result from [11] 
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The thermal resistance for a stripe is from [12], see Figure 4. 
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In regions where the thermal fluxes intersect, the formulas from [13] (with typos 
removed) are applicable. 
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The step-by-step procedure starts with the heat sources on the top surface; which 
are strips with the short side dimension being that of the gate, or the gate-drain spacing. 
Using eqns. (8) and (9), to determine the flow pattern; the footprint on the next lower 
surface is determined. The thermal resistances may be obtained from Eq.(10) which is 
applicable for either a single strip, or multiple ones; whereas eq.(7) is restricted to single 
ones. With the flow patterns and resistances of each layer determined, the calculation 
continues by starting at the assumed sink temperature and working upward through the 
The change in temperature across a layer is given by 

- the temperature rise 
= the thermal resistance 
= the power dissipated 
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The actual temperatures are used in the layers with constant thermal conductivity, and the 
artificial ones are used in the GaN, SiC, and Sapphire layers. For the GaN, SiC, and 
Sapphire layers, the actual temperature is determined from the transformed one by 

ffrf - -9.x 
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While the KirchhofF transform gives a linear equation, the nonlinearity inherent in eq.(l) 
is now pushed into the boundary conditions. For a layered structure such as a HEMT, 
this can be troublesome when attempting to write simple code [14]. In layered structures, 
the artificial temperature will generally have jumps across interfaces. However, as stated 
in [14], these jumps can be eliminated by choosing the analytical forms for the 
conductivities of the layers, (f), in such a way that their ratio is independent of 
temperature. This lead us to assume the forms given above. For GaN and SiC, the ratio 
of the conductivities is 
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or a 6% change over a 200 C range. 
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III. Comparison with Experimental Data 

Temperature measurements on AlGaN/GaN HEMTs on both SiC or Sapphire 
were made using Raman spectroscopy [15], [16]. For a single gate device [15] on SiC, 
the measured temperature on the top surface was 1 80 °C. We choose the heat source to 
be 1/i. wide (the gate width), and the GaN thickness was F - \ 2p- . Using eq.(10) for n 
= 1 gives Ql: ft = 0.587. Assume X (t) =1.6 (the 300 K value) as we don't 
know the actual temperature at the interface yet. This will underestimate the final 
channel temperature. We find 0(6oW) - /8.2S K /w . Using eqns.(7) & (8) we find the 
footprint is 3.086 by 200.2 jll . Repeat eq.(10) and obtain 6^ ft - 2.02, and 
0(i/c) = 3.7- %/ . The dissipated power is 1.754 W, which yields AT(sjc)= 52 K, 

AT = 32.2 K. Using eq.(15) gives T(siC y adfaJ.) - 86.6 °C. Then in the GaN, T e 
= 359.6 + 32.2 = 391.8 K, or T(g*h, oaJmM. ) - 121 °C. The measured value was 

124°C - 5 °C. We can make a correction for by using X(359.6 K) = 1.24. Then 

^(Ga N) - 23.7 KAV and the temperature becomes 132° C. The 3-D model used in [15] 
assumed constant ft and predicted T = 99° C. For the sapphire substrate case we choose 
L = 4fL (gate length), 3ook) = 0.49; then Q(^^j) - 8.35 K/W, = 180.5 

K/W, A T(Sofy) - 1 17.3 K, = 5 4 Then the temperature on the top surface 

becomes 176°C. The measured value was 180°C; and the 3-D model predicted 140 °C. 

An 8-gate device on SiC (for 2-4 GHz amplifier applications) was measured in 
[16], and the thermal profile is shown in Figure 5 a). Using the outlined procedure one 
obtains T = 179 C (shown dotted); which agrees favorably with the average of the 
measurement. 

Measurements using nematic liquid crystal thermography on a 2-gate device on 
sapphire were performed in [17]; see Figure 5 c). Our computations yielded 59 C (lower 
curve) and 76.5 °C for the upper case. 



IV. Conclusions 

We have demonstrated an easy method to estimate the channel temperature in an 
AlGaN/GaN HEMT for microwave amplifiers. With accurate expressions for the 
temperature dependent thermal conductivity and using previously published equations for 
the thermal resistance of parallel stripes, one calculates a temperature that is within 5- 
10% of the true average channel temperature. 
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Figure Captions 

1 . Schematic of an 8-gate HEMT for Ka-band operation. 

2. Approximate layers of material with average, or ranges, of thermal conductivities. 

3. Device of Figure 1 on an approximate MMIC chip. The thermal resistance of 
each layer is calculated separately (9 = #r A )• 

4. Geometry for calculating thermal resistance for a single stripe a), or parallel 
stripes b). 

5. Comparison with Kuball [16] and Park [17]. Part a) shows measurement (solid) 
[16] and estimate (dotted). Part b) shows estimate (dotted) for two power 
dissipation cases [17]. 
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